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ABSTRACT 

A  compensation  scheme  is  presented  for  the  pitch-yaw 
control  system  of  a  strategic  missile  which  maintains 
stability      and    performance      while  the      missile    spins.  The 

conventional  missile  control  design  for  boost  has  uncoupled 
pitch,  yaw  and  roll  channels  ideally  suited  for  a  non- 
rolling  missile.  The  effect  roll  has  on  the  conventional 
pitch-yaw  controller  is  discussed  and  the  development  and 
simulation   of    a   modified   system   is   presented. 
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I.    IHTIQDOCTIQH 

Future  ballistic  missiles  nay  need  the  ability  to  'spin 
during  boost  ascent  through  the  atmosphere  as  a  counter- 
measure  against  first-generation  continuous-duty  laser 
threats.  Since  the  conventional  missile  control  design  has 
uncoupled  pitch,  yaw  and  roll  channels,  some  pertormance  and 
even  system  stability  may  be  lost  in  rolling  the  missile. 
This  study  evaluates  the  performance  of  a  conventional 
pitch-yaw  controller  in  the  presence  of  spin.  A  compensa- 
tion    scheme   is   then    proposed   to   restore   system    performance. 

The  block  diagram  of  a  submarine- launched  ballistic 
missile's  pitch-yaw  control  system  is  used  as  the  design 
vehicle  for  this  study.  Since  the  physical  parameters  of 
this  system  are  classified,  only  the  block  diagram  is  used 
in  the  analysis  and  design  process.  All  numerical  values 
needed  tor  analysis  and  simulation  are  borrowed  from  other 
systems    or  derived   as   needed. 

The  study  is  conducted  as  follows.  First,  the  basic 
laws  of  motion  and  thrust  vector  control  are  introduced. 
Together,  these  sets  of  eguations  define  the  flight  dynamics 
of  the  missile  and  the  method  for  controlling  the  missile. 
The  pitch- yaw  controller  for  the  non-rolling  missile  is  then 
designed  and  simulated.  This  system  has  uncoupled  channels 
for  pitch  and  yaw  control  which  become  coupled  when  the 
missile  is  rolled.  Next,  the  missile  is  forced  to  roll  and 
the  effect  cross-coupling  has  on  the  system*s  performance  is 
analyzed.  A  compensation  scheme  is  then  proposed  to  restore 
the  system's  perf ormance.  Analysis  and  design  are  carried 
out  using  classical  and  modern  controls  analysis  and 
simulation   results  are    presented. 


A.   HISSILE  D1HAH1CS  AMD  THBUST  VECTOR  CCMIBOL 
1.   Notation  and  Reference  Axis  JLystei 

The  reference  axis  system  used  tc  define  vector 
quantities,  dimensions  and  angles  of  interest  is  shown  in 
Fig-  1-1-  The  axis  system  is  fixed  within  the  missile  body 
with  its  origin  at  the  missile's  center  of  gravity.  The 
orthogonal  axes,  x,  v  and  z,  will  also  be  referred  to  as  the 
roll,  pitch  and  yaw  axes,  respectively.  Pig.  1-1  indicates 
that  the  right-hand  rule  is  used  to  determine  the  positive 
direction  for  rotational  vector  quantities. 
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Fig.  1-1   R€ference  Axis  System 


Table  1-1  lists  the  symbols  used  to  refer  to  key 
parameters  and  quantities  of  interest.  Other  symbols  «ill 
be  defined  as  they  arise.  lhe  moments  of  inertia  about  the 
center  of  gravity  (eg)  are  defined  as: 


A  =  2  6m    (y2*z2) 


(1.1) 


B  -  2  da  (z*+x*) 


(1.2) 


C  =  2  d"m  (x**y*) 


(1.3) 


The  products  of  inertia  are  defined  as 


D  =  2  (<fm)  yz 


(1.4) 


E  =  2,  (6m)    xz 


(1-5) 


F  =  £  (dm)  xy 


(1.6) 


Where  £6m   symbolizes  the  sum   of  differential  mass  elements 
which  make  up  the  missile  tcdy. 


TABLE  1-1 
NOTATION  EEFINITIONS 


1 

Roll  Axis 

X 

Pitch  Axis 

y 

, , 

Yaw  Axis 

z 

Angle  of  incidence 

X 

0 

V 

Angular  rate 

p 

q 

r 

Component  of  missile 
velocity  along  axis 

U 

V 

w 

Component  of  force 
on  missile  along  axis 

X 

Y 

Z 

Moments  acting  on 
missile  about  axis 

L 

M 

N 

Moments  of  inertia 
about  each  axis 

A 

B 

C 

Products  of  inertia 

D 

E 

F 

..  .. 

2-   Euler fs  Equations  of  Motion 

The  equations  of  notion  foe   a  body  with  six  degrees 
of  freedom  and  a  mass,  id,   are  summarized  below,        A  dot  (•) 
over  any  symbol  refers  to  the  time  rate  of  change,  d/dt,   of 
the  guantity  indicated-   [fief.  1] 


m(U  +  qv    -    rv)  =  X 


(1.7) 


a(v  ♦  rU  -  pw)  =  1 


(1.6) 


n (w  -  gU  +  pv)  =  Z 


d-9) 


Ap  -  (B-C)qr  ♦  C(r2-g2)  -  E(pq*r)  ♦  f  (rp-q)  =  L    (1.10) 


Bq  -  (C-A)rp  ♦  E  (p2-r2)  -  F(gr*p)  ♦  D  (pg-r)  -  H    (1-11) 


Cr  -  (A-B)pq  +  F  (q2-p2)  -  D(rp  +  q)  ♦  E  (qr-p)  =  N    (1.12) 


Ihe  first  three  equations  above  are  force  equations 
and  the  second  three  are  moment  equations.  In  particular, 
Eqs.  (1-11)  aod  (1.12)  are  tbose  from  which  pitch  and  yaw 
angular  accelerations  are  found-  Note  that  the  second  term 
in  each  of  these  equations  constitutes  a  cross-coupling  term 
between  the  two  channels  which  is  proportional  to  th€  roll 
rate.  This  cross-coupling  effect  is  eliminated  in  the 
conventional  pitch-yaw  controller  design  by  preventing  roll 
and  thus  setting  the  roll  rate,  p,  equal  tc  zero- 
Assuming  that  the  missile  body  is  approximately 
rigid  leads  to  the  conclusion  that  the  products  of  inertia 
are  equal  to  zero.  Thus  Egs.  (1.11)  and  (1.12)  are  reduced 
tc: 


Bq  -  (C-A)rp  =  M  (1.  13) 


Cr  -  (A-B) pq  =  N  (1-  14) 


3-   Thrust  Vector  Control 

Ihe  only  means  used  to  control  the  flight  of  the 
missile  under  study  is  to  deflect  the  exhaust  nozzle  as 
shown  in  two  dimensions  in  Fig.  1-2.  Deflecting  the  exhaust 
nozzle  causes  the  driving  force  of  the  rocket  motor  tc  fce 
pointed  in  a  direction  other  than  along  the  missile's  longi- 
tudinal axis.  Ihis  action  creates  force  components 
perpendicular  to  the  roll  axis  which  do  not  act  through  the 
missile's  center  of  gravity.  The  resulting  moments  are 
about  the  y-  and  z-axes  and  cause  the  oissile  to  pitch  and 
yaw.  The  moments  about  the  center  of  gravity  caused  ty  the 
deflected  thrust  vector  are: 


H  =  <TZ)1C  (1-15) 


N  =  {1/)lc  (1.16) 


Where  lc  is  the  length  from  the  missile's  center  of  gravity 
to  the  pcint  of  application  of  the  thrust  (exhaust  nozzle) 
and  1y  and  lz  are  the  lateral  thrust  components  parallel  to 
the  y-axis  and  z-axis.  If  dg  and  dy  are  the  angles  of 
nozzle  deflection  in  the  pitch  and  yaw  planes,  then  the 
perpendicular  thrust  components  are  related  to  the  tctal 
rocket  motor  thrust  by  the  following  eguaticns: 


ly    =  I-sin(dy)  (1.  17) 


1Z    =  T«sin(d0)  (1.18) 
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x-axis 


Fig.  1-2   Thrust  Vector  Control 

Assuming  that  the  maximum  angle  of  nozzle  deflection 
is  snail,  as  it  must  be  for  any  physical  system  of  this 
type,  permits  further  simplification  of  Egs.  (1.15)  and 
(1.16)  to  the  lin€ar  Egs.  (1.21)  and  (1.22). 


sin(de)  ^  de 


(1.19) 


sin  (duy)  =  dy 


(1.20) 
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a  =  i-a0-ic 


(1.21) 


N  =  l'dy-lc 


(1.22) 


Equations  (1.21)  and  (1.22)  give  the  moaents  in 
pitch  and  yaw  caused  by  the  cocket  motor.  Combining  these 
two  equations  with  the  general  rotational  motion  eguations 
in  pitch  and  yaw,  Eqs.  (1.13)  and  <1.14),  permits  formula- 
tion of  the  general  differential  equations  for  pitch  and 
yaw. 


T-de-lc  =  Bg  -  (C-A)rp 


(1.23) 


T*dv|/*lc  =  Cr  -  (A-B)pg 


(1.24) 


Solving  for  the  angular  acceleration  components  gives: 


q  =  (1/B)£l-de-lc  ♦  (C-AJrc] 


(1-25) 


r  =  (1/C)[T-df  -lc     ♦  (A-E)pq] 


(1.26) 


Equations  (1.25)  and  (1.26)  point  out  the  cross- 
coupling  effect  that  the  roll  rate,  p,  has  on  the  system. 
For  instance,  if  the  roll  rate  is  set  equal  to  zero  the 
equations  for  pitch  and  yaw  are  completely  uncoupled. 


q  =  (1/E)T«d9-lc 


(1-27) 


r  =  (1/C)  T«dy  lc 


(1.28) 
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Figure  1-3  depicts  the  rolling  pitch-yaw  system  just 
described  in  block  diagram  form  and  highlights  the  cross- 
coupling  terms  introduced  by  roll.  Ihe  classical  missile 
control  system  design  approach  assumed  that  the  roll  rate 
was  small  and  went  to  great  lengths  to  ensure  this  was  true. 
Making  the  assumption  that  the  roll  rate  is  negligible 
allows  the  system  designer  to  separate  the  pitch  and  yaw 
channels  and  deal  with  them  independently.  This  study  will 
analyze  the  effect  cross-coupling  has  en  stability  and 
performance  and  formulate  a  compensation  scheme  for  the 
cross-coupled  system. 
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II.    fiEVELOPgENT    OP    TgB    NOM -BOLL IMG    HISSILE 

The  strategic  Bissile  whose  pitcn-yaw  control  systea 
will  be  used  as  the  design  vehicle  for  this  study.  The 
block  diagram  of  the  system's  pitch-yaw  control  system 
(without  roll)  is  shown  in  Fig.  2-1.  Since  the  physical 
parameters  of  this  missile  are  classified,  only  the  block 
diagram  will  be  used.  All  system  constants  and  physical 
parameter  values  which  are  used  for  analysis  and  simulation 
will    be   taken   from  other  systems   or  derived   as    needed. 

The  missile  is  a  three-stage  booster  rocket.  Flight  of 
the  missile  is  controlled  by  vectoring  the  exhaust  thrust  as 
described  in  the  preceding  section.  Based  upon  predeter- 
mined flight  trajectories  and  measured  or  estimated  parame- 
ters during  flight,  the  guidance  computer  generates  pitch 
and  yaw  commands  and  transmits  these  to  the  pitch-yaw 
controller      [fief.    2].  The   algorithm      which      the      guidance 

computer  uses  to  generate  the  pitch  and  yaw  commands  will 
not      be   addressed.  Only      the      response    of      the      pitch-yaw 

controller    to   a    general   set    of   commands    will   be   considered. 

Table  2-1  lists  the  physical  parameter  values  used  at 
the  various  stages  for  design  and  siuulation.  These  parame- 
ters tear  no  reseablance  to  those  of  the  Trident  II  missile 
but  serve  only  to  provide  actual  numbers  for  the  variables 
so  that  simulation  could  be  done.  The  paraaaters  listed  in 
Table  2-1  were  actually  taken  froa  the  Saturn  V  Rocket;  also 
a  three-stage  booster  rocket.  The  moments  of  inertia  listed 
in  Table  2-1  were  obtained  from  the  other  parameters  by 
assuming  that  each  stage  has  approximately  the  same  inertial 
characteristics  as  a    uniform   right-circular   cylinder. 
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Fig.    2-1      Non-Kolling   Pitch-Yaw   Control    System 

Figure  2-1  shows  the  schematic  block  diagram  of  the 
system  used  to  control  pitch  and  yaw  in  the  missile.  Since 
the  missile  does  not  spin,  the  pitch  and  yaw  channels  are 
uncoupled.  The  system  consists  of  a  servo- act uator  used  to 
deflect  the  exhaust  nozzle,  the  missile  dynamics  and  sensors 
which  feed  back  angular  velocity  and  position.  The  deflec- 
table exhaust  nozzle  saturates  at  a  maximum  deflection  of  ♦ 
five  degrees.  In  this  analog  system  all  gains  are  constant 
throughout    operation. 
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TABLE  2-1 


PHYSICAL  PABAflETEE  SPECIFICATIONS 


1st  Stage        2nd  Stage          3rd   Stage 

Initial   Mass                     slugs 

1.488x105 

3.22x104 

8.14x105 

Mass  Rate                           slug/sec 

930 

81.49 

14.75 

Thrust                                  lb- 

7.496x106 

1.125x10 

2.249x105 

Length                                  ft 

135 

130 

100 

Burn  Time                          sec 

150 

359 

479 

■ 

Moment  of    Inertia            ,         _  2 
about  Pitch  Axis       sLnS-it 

m(1542.1J 

m(l476.4) 

m(858.3; 

Moment  of    Inertia         n         „.2 
,       .   v        .     .          slug-ft 
about  Yaw  Axis 

m(l542.l) 

m(l476.4; 

m(858.3) 

Moment  of    Inertia            ..        r+2 
about  Roll  Axis          s^g-ix 

m(i36.i; 

m(136.i; 

m(50.0; 
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A.        DESIGN    OP    THE    SEfiVO-ACT OATOB 

The  transfer  function  of  the  servo-ac tua tor  which 
deflects  the  exhaust  nozzle  of  the  missile  is  given  in  Eg. 
(2.1).  The    task      is   to      select   nuaerical      values   for      the 

parameters,  a,  b  and  f,  such  that  the  servo-actuator  has 
acceptable  performance  characteristics.  Again,  the  param- 
eter values  will  be  selected  without  regard  to  those  used  on 
the    actual    missile. 


dp  (s)      _.       abz         (2.1) 

E,  (s)  (s   ♦    a)  (s*   ♦    2f  bs    ♦   b^) 

It  is  desired  that  the  actuator's  performance  resemble 
that  cf  a  well  damped  second-order  system  with  linimum  over- 
shoot and  a  short  settling  time  in  response  to  step  inputs. 
These  performance  characteristics  need  to  be  more  stringent 
than  those  imposed  on  the  entire  system  so  that  the  actuttor 
does  not  become  a  limiting  element  when  the  entire  system  is 
considered    later.  The  design   specifications    used      for    the 

actuator  are   summarized   as    follows. 

1.  Pole  Locations 

In  order  to  ensure  that  this  third-order  system 
behaves  essentially  as  a  second-order  system,  the  real  pole 
must  have  a  magnitude  significantly  greater  than  the  real 
part    of    the   complex    roots  of   the  second-order   term,    e.g. 

a      >    10fb  (2.2) 

2 .  Damping 

In  order  to  make  the  dominant  complei  pair  of  poles 
exhibit    good    damping,    the    value   of    £       must   approach    one. 

1    >  f    >    0.8  (2.3) 
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3.   Settling  lime 

In  order  to  keep  the  actuator  from  Halting  the 
settling  time  of  the  entire  system,  the  servos  settling 
time  must  be  significantly  shorter  than  that  of  the  pitch- 
yaw  controller-  If  the  desired  settling  time  of  the  entire 
system  is  to  be  approximately  ten  seconds  and  the  actuator's 
settling  time  is  3%  of  the  system's,  then  the  actuator's 
settling  time  will  be  selected  according  to  the  relationship 
below. 


t      .,,  ,  <    0.03    t      ,  .,  .  .  (2.4) 

settling  servo  -  settling  system  *         ' 

4.      leak    Cvershcot 

Ihe  peak  overshoot  of  the  servo-actuator  must  be 
small.  Otherwise,  it  may  have  an  appreciable  effect  on  the 
overall   system's    peak    overshoct- 


M(servc)    -    1    <    1% 


(2.5) 


Where  a  (servo)  is  the  peak  overshoot  of  the  actuator. 

Equation  (2.5)  permits  selection  of  the  appropriate 
damping  ratio,  f.  If  the  system  is  assumed  to  exhibit 
second-order  characteristics,  then  Eg.  (2.6)  can  be  solved 
for  the  damping  ratio. 


M  (servo)  -  1  =  exp 


-  1T£_ 
vT^"f 


<  0.01 


(2.6) 


Solving  for  the  damping  ratio,  f,  yields  the  following. 


f  > 


[ln(0.01)  ]* 


=  0.826 


r2  ♦  [lu(O.OI)  ]2 


So,   a  damping   ratio  of  0.9  Kill  be   used. 
also  satisfies  Eg.  i2»3). 


(2.7) 


This  selection 
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The  specification  of  Eg.  (2.1),  combined  with  the 
above  choice  of  damping  ratio,  can  be  used  to  select  the 
natural  freguency  of  the  second-order  tecs. 


settling  servo 


4_ 
fb 


<  0.3  sec 


(2.8) 


Solving  for  the  natural  freguency  leads   to  a  selection  of  b 
egual  to  15. 


b  >  14.8  rad/sec 


(2.9) 


Equation  (2.2),  combined  with  the  above  choices  of 
damping  ratio  and  natural  freguency,  leads  to  a  selection  of 
the  real  pole,  a,  egual  to  150. 


a  >  10ft  -  135 


(2-10) 


Ihe  root  locus  of  the  system^  characteristic  egua- 
tion,  obtained  by  varying  the  real  pole,  a,  is  shown  in  Fig. 
2-2.  Highlighted  are   the  root   locations   for  the   above 

choice  of  parameter  values.  As  shown  in  Pig.  2-2,  ail  the 
roots  are  well  into  the  left  half  of  the  s-plane,  a  large 
phase  margin  exists  for  the  second-order  conplex  pole  pair 
and  the  second-order  pair  is  far  enough  to  the  right  cf  the 
single  real  pole  to  achieve  dominance.  Figure  2-3  shows  the 
open-loop  freguency  response  of  the  actuator  and  lakes  the 
large  phase  margin  more  evident.  Figure  2-4,  the  closed- 
loop  freguency  response,  points  out  the  linitation  of  the 
actuator  to  follow  input  signals  of  freguency  greater  than 
about  10  radians/second.  Figure  2-5  shows  the  step  response 
of  the  system  and  as  indicated,  it  meets  the  design  specifi- 
cations. Note  that  although  the  simulation  shown  in  Fig. 
2-5  is  that  cf  a  full  step  of  unit  magnitude,  the  nozzle 
actuator  saturates  at  a  maximum  angle  of  deflection  of 
0.0873  radians  (5  degrees). 
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Fig.  2-2   Servo-Actuator  Boot  Locus  Plot 

Figure  2-2  shows  that  the  magnitude  of  the  real  pole 
is  far  greater  than  the  real  part  of  the  complex  pair  and, 
therefore,  will  have  little  effect  on  the  system's  perform- 
ance. Also,  the  dominant  seccnd-order  pair  is  well  damped 
with  approximately  a  72-degree  phase  margin  and  an  effective 
damping  ratio  of  0.95.  Ihe  system  is  static  as  long  as  the 
term,  ab2,  remains  less  than  7.57  x  105. 
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B.   DESIGN  OF  THE  PITCH-XAH  COKIBOLLER 

The  block  diagran  of  the  pitch-yaw  control  systea  is 
shown  in  Fig.  2-1-  Since  the  pitch  and  yaw  channels  are 
identical  and  decoupled  only  one  of  them  need  be  considered. 
The  paraaeters  K  (K  )  and  K  (F+)  will  fce  selected  so 
that  the  systea  will  exhibit  ainimal  overshoot  to  step 
inputs  and  have  a  settling  tiae  less  than  ten  seconds.  The 
one-channel  systea  to  be  considered  is  shown  in  Fig.  2-6. 
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Fig.  2-6   Fitch  or  law  Channel  Block  Biagrai 


In  order  to  begin  the  analysis,  it  is  necessary  to 
temporarily  ignore  the  nonlinear  saturation  of  the  servo- 
actuator.  The  linearized  transfer  function  for  the  pitch- 
yaw  control  systei  then  becomes  that  of  Eg.  (2.11). 

e(s)      =  Keab2Tlc/E  (2.11) 

9c(s)  ss    ♦   s*(a*2fli)     ♦   s3(b*«-2fba)     +    s2  (ab2) 

♦    s(KreKQafc2Tlc/E)    ♦    Keab211c/E 

The  denominator  of  Eg.  (2.11),  when  set  egual  to  zero, 
defines  the  systeu*s  characteristic  equation.  The  charac- 
teristic eguation  nay  be  partitioned  as  follows.  Define  a 
new  variable,  J. 


J  =  ab2K(Tl  /B)  (2.  12) 

e   c 

The  characteristic  eguation  now  becomes; 


ss  ♦  s*(a  +  2fb)  ♦  s3(b*+2ffca)  (2.13) 

♦  s2  (ab2)  ♦  s  (Ftq  J)  +  J  =  0 

Hearranging  yields  the  partitioned  form  shown  below. 

Kres  +  1  =  -1        (2.  14) 

s»  ♦  s«  (a  +  2f b)  ♦  s3(b2  +  2fba)  ♦  s2(ab2)     J 

Eguation  (2.14)  is  in  the  form  for  root  locus  analysis 
where  J  is  treated  as  the  open-loop  gain  and  the  left  side 
of  the  eguality  as  the  open-loop  transfer  function.  Varying 
Krs  and  J  together  sweeps  out  a  region  in  the  s-plane  into 
which  the  roots  of  the  characteristic  eguation  can  be 
located.  These  "root  relocation  zones"  are  shown  in  fig. 
2-7.   [Bef.  3] 

Figure  2-7  was  made  by  varying  J  and  K re  in  Eg.  (2.  14) 
and  plotting  the  root  locations.  For  example,  if  Kre  =  0 
then  Eg.  (2.  14)  contains  five  finite  poles  and  five  zeros  at 
infinity.    The  root  locus  will  follow  the  path  indicated  in 
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Pig.  2-7  (as  J  increases).  Note  that  the  real  pole  at  -150 
is  not  shown  in  Fig.  2-7.  As  Kre  is  increased  incrementally 
from  zero  towards  infinity  and  the  sate  procedure  is  used, 
i.e.,  plot  the  root  locus  over  a  large  (positive)  range  of 
J,  then  a  family  of  root  loci  will  be  swept  out.  The  cross- 
hatched  area  in  Fig.  2-7  is  the  area  enlosed  by  the  root 
loci  for  J  and  KTO  varying  from  0  to  co. 


-60.0    -60.0     -40.0     -30.0     -80.0     -10.0 


0.0 


10.0      20.0 


Fig.  2-7   Boot  Relocation  Zones  of  Pitch-Xaw  Controller 
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Figure  2-7  indicates  that  the  system  is  unstable,  with 
two  roots  in  the  right  half-plane,  for  values  of  Kre  less 
than  or  egual  to  zero-  The  values  of  K  re  and  J  are  esti- 
mated as  0.8  and  1.57  X  105,  respectively.  This  selection 
locates  the  characteristic  eguation  roots  as  shown  in  fig. 
2-8.  Ihe  open-loop  gain,  K0  ,  was  found,  from  the  definition 
of  J,  to  range  between  the  values  shown  below. 

1st  Stage  2nd  Stage  3rd  Stage 

1.071-0.07  1.51-0.  14  1.45-0.19 

Exact  solutions  for  Ke  and  K  re  were  not  obtained  analyt- 
ically since  the  system  contains  a  nonlinear  saturation 
effect  which  was  ignored  in  constructing  the  root  locus. 
Also,  Ke  is  a  function  of  system  parameters  which  change 
continuously  during  flight  and  dramatically  at  booster  stage 
separations.  To  fine  tune  the  selection  of  Ke  and  Kre  , 
simulation  studies  were  conducted  using  DSL/360  during  each 
of  the  three  booster  stages  of  operation.  Ihe  results  of 
the  simulations,  for  various  values  of  F  ,  are  shown  in 
Figs.  2-9,  2-10  and  2-11,  where  8  is  given  as  a  function 
time. 

The  simulations  illustrated  that  the  gain,  f^  ,  has 
little  effect  on  the  system's  performance.  This  is  because 
it  is  "upstream"  from  the  saturated  nozzle  actuator  and,  as 
such,  can  only  shorten  the  tine  required  for  the  actuator  to 
reach  maximum  deflection.  Since  this  time  is  already  small 
compared  to  the  system's  response,  Ke  has  negligible  effect. 
To  simplify  the  problem,  K9  was  set  equal  to  one.  The 
velocity  feedback  gain  Kre  has  a  pronounced  effect  on  both 
the  overshoot  and  the  settling  time  of  the  system.  from  the 
simulations,  Kre  was  selected  to  be  1.5.  This  choice 
ensures  the  specified  settling  time  of  ten  seconds  during 
all  stages  and  results  in  overshoot  only  during  stages  one 
and  two,  where  it  is  less  than  ten  percent. 
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Fig-  2-8   Rcot  Selection  of  Partitioned  Equation 
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TIME  (SEC) 


Fig.  2-9   First  Stage  Simulation  of  Pitch-¥a¥  Controller 
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Fig.  2-10   Second  Stage  Siaulation  of  Pitch-law  Controller 
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Fig.  2-11  Third  Stage  Simulation  of  fitch-Yaw  Controller 
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III.  gFFBCI  OP  BOLL 

Rhen  the  non-rolling  pitch-jaw  control  system  just 
designed  is  forced  to  roll,  the  system  becomes  cross-coupled 
as  shown  in  Fig.  3-1.  The  mechanism  which  creates  and  main- 
tains the  rolling  motion  will  not  be  considered.  It  is 
assumed  that  a  control  system  outside  the  pitch-yaw 
controller  and  independent  of  the  main  thrust  motor  main- 
tains an  approximately  constant  roll  rate,  p.  Furthermore, 
it  is  assumed  that  the  roll  rate  can  be  measured  during 
flight  or  estimated  from  treasured  parameters  with  reasonable 
accuracy.  Ihe  purpose  of  this  section  is  to  analyze  the 
effect  roll  has  on  the  previously  non-rolling  pitch-yaw 
controller. 

The  Laplace  transfornations  of  the  linearized  non- 
rolling  system's  governing  differential  eguations  can  be 
written  as  two  separate  eguations;  one  for  the  pitch  channel 
and  one  for  the  yaw  channel. 


s*9  ♦  GK  res6  ♦  G6  =  G9C 


(3.1) 


s^y  ♦  HK^sf  ♦  H?  =  HY, 


lrv 


13.2) 


where  G  and  U  are  defined  as 


G    =    RJTL/B) 


ab* 


L   (s*a)  (s*  +  2fts*b2)-l 


{3.3) 


H    =    Ky<Tlc/C) 


ab* 


L  (s*a)  (s2*2fbs  +  b2)-l 


{3.14) 


32 


q> 


+ 

co 
I 


• 

t3 

M 

o 

+> 

z$ 

N  +» 

N    O 

o  < 

SZ     1 

o 

t 

<X1 

CO 

i 

1 

o 

M 

CO 


CO 


•  a* 


+  T  + 


o 

El  I  w 


O 


i 


i 

o 


CQ 


PP 
I 


O 


CO 


t-|CO 


•  M 


o 


o 


^ 

Nozzle 
Servo -Actuator 

i 

1 

3- 

A 


v<J 


CO 


M 
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Equations  (3.1)  and  (3.2)  can  be  rearranged  and  put  into 
matrix  form  as  shown  in  Eg.  (3.5).  Note  that  the  eguation 
has  only  diagonal  non-zero  elements  and  isr  therefore, 
uncoupled. 


"s2*GKres*G       0 

0       s2*HKr*s  +  HJ 


9 

- 

G 

0 

9c 

r 

.0 

H. 

_Tc. 

(3.5) 


Rolling  the  missile  leads  to  the  addition  of  cross- 
coupling  terms  between  the  equations  governing  pitch  and  yaw 
response.  The  transformed  differential  eguations  governing 
the  responses  of  the  rolling  system  are: 


s*0   ♦  GK^  s6  +  G8  +  gsV  -    G0C 


(3.6) 


S*T      ♦    HKrysV    ♦    HY    ♦    Bs8    =    Efc 


(3.7) 


Where  G  and  H  are  defined  in  Egs.   (3,3)    and  (3.<4)  and  C  and 
R  are  defined  in  Egs.  {3.8)  and  (3.9). 


C  =  (A  -  C)p/E 


13.8) 


R  =  (E  -  A)p/C 


(3.9) 


Rearranging  Egs.   (3.6)  and  (3.7)  into  matrix  form  as  before 
leads  to  the  coupled  form: 


32  +  GKj.q  s*G 
£s 


s2*HKrfs*H 


0 

= 

G       0^ 
0       Hj 

(3.10) 


Eguation  (3.10)  shows  that  the  introduction  of  roll  has 
affected  the  transfer  functions  which  determine  the  input- 
to-output  relationships  of  the  system   in  twc  ways.    first, 
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it   has   changed   the   direct  path   transfer   functions   of 

pitch-to-pitch  and   yaw-to-yaw  command-to-cutput.     Second, 
roil  has  introduced   pitch-to-yaw  and   yaw-tc-pitch  comand- 

to-output   transfer  functions  which   sere   zero  before   the 

missile  was  forced  to  roll.  Ihe  effect  that  these  changes 

have  on  the  systea's  stability  and  performance  will  be  exam- 
ined next, 

A.   EFFECT  OF  BOIL  OH  SYSTEH  STABILITY 

As  pointed  out,  the  rolling  system  is  governed  by  four 
transfer  functions  (2  inputs  x  2  outputs) .  Consider  cnly 
the  pitch-to-pitch  input-to-output  transfer  function,  i.e., 
hold  the  yaw-command  input  equal  to  zero  and  ignore  the  yaw 
output-  Under  these  conditions,  the  block  diagram  of  the 
system  may  be  rearranged  as  shown  in  Fig.  3-2.  The  transfer 
function  obtained  from  the  diagram  is  given  in  Eg.  (3.11). 
As  indicated,  the  pitch  response  of  the  system  is  influenced 
by  a  variety  of  system  parameters,  as  it  was  before  the 
missile  was  rolled.  Now,  however,  the  root  locations  of  the 
system  are  also  affected  by  the  roll  rate. 

8(S)  ■  N5s5+N4S**K5s3+N2s2*N1s+N0  (3.11) 

9c(s)       s»o*Cgs9*D8s8+DyS7*D6s*+DcSS4D4s4 

♦  D,  s3«-E^  s2  +  L,  s  +  Dq 

where ; 

N5    =     (abVB)TlcKe 


UA    =    (ab2/B)llcKe(a*2ffc) 


N5    =  (ab2/B)TlcKe(bs  +  2fba) 

N2    =  (ab2/B)IlcK9(ab*) 

N1    =  (ab2/B)TlcF9(FrvFvab2Tlc/C) 

Nn    =  (ab2/B)Tl    K     (Ku/abZTL, /C) 
u  c    e       >  c 

D^     =  2  (a  +  2fb) 
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D   =    2(b2*2fba)     ♦    £  (A-B)  (A-C)  p2/CB  J    ♦    (a*2ft)2    =    D81     ♦    D82  p* 
Dy=    2ab2    ♦    2(a*2fb)  (b2  +  2fta  +  [ <A-B)  (A-C)p2/CB])     =    D71     ♦    C72  p2 

D6=    (b2*2fba)  (b2*2fba    «•    2[  (A-E)  (A-C) p2/CE J) 

♦  (ab2Tlc)[  (KryKt/C)+|^.e  K9/B)  J    ♦     (a*2ab)  (2ab2) 

♦  [  (A-B)  <A-C)p2/CB](a*2U)2    =    D61      ♦    %2   p2 

D5=    (ab2Tlc)[  (Ke/BJ  +  lKv/CJ  ]    ♦    2[  ( A-B)  (A-C) p2/CB  ]  (ab2) 

♦  2(b2+2fba)  (ab2)     ♦    2[ (A-B)  (A-C) p2/CE  ]  (b2  +  2f ba)  (a*2f b) 

♦  (a»2fb)  (abSTJcH  (KrH,  Kv/C)  ♦  (Kr9  Ke/E)  ]    =    D51    4    C52  P2 

B4=    (at2Tlc)[  (Kv/C)  +  (Ke/E)  ](a+2fb) 

+    2[  (A-B)  (A-C) p2/CB]  (a*2f b)  (ab2) 

♦  [  (KryKr/C)»  (Kre  Ke/E).]  (ab2Tlc)  (b2*2fba) 

♦  [  (A-B)  (A-C) p2/CB](b2*2f ba)2    ♦    (at2)2   =    D41    ♦    G42P2 

D  5  =    [  (Ky/C)  ♦  (Ke/B)  ]  (abzi^  )  (t2+2fba) 

♦  2£  (A-B)  (A-C) p2/CB]  (at2)  (b2*2fba) 

♦  [  (KrrK*/C)     ♦    (^  Ke/E)  ]((ab2lJc  )at2)    =    D  31    ♦    13  32P2 

D2     =    [  (Kw/C)*(Ke/E)  ](ab2)  (ab2llc) 

♦  [  (A-B)  (A-C)p2/CE](at2)2 

♦  (K^K^K^  Ke/BC)  (Tlcat2)2    =    D21     ♦    D22  p2 

D-!     =    Ke  (ab2Tlc/CB)  2Kr+,  Kt    ♦    Kre  K  e(ab2Tlc/CE)  zr^ 
D0    =   £  {!W  <ab2llc)2]/CB 
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The  denominator  of  Eg-  (3.11),  when  set  equal  to  zero, 
defines  the  characteristic  equation  of  the  system. 
Examination  of  Eg.  (3.12)  (in  which  Eg.  (3.11)  defines  the 
coefficients)  shews  that  as  the  roll  rate  is  increased  from 
zero,  the  roots  of  the  characteristic  equation  will  move  as 
a  function  of  the  sguared  rcll  rate.  For  analysis,  the 
region  of  operation  will  re  restricted  to  the  beginning  and 
end  of  each  booster  stage.  At  these  instants,  the  charac- 
teristic eguation  is  as  given  in  Eg.  (3.12)  with  numerical 
values  for  th€  coefficients  given  in  Table  3-1.  Ihe  end  of 
stage  three  is  not  considered  because  the  missile  is  no 
longer  performing  powered  flight  at  the  end  of  the  final 
stage. 


s»o  ♦  D9S»  ♦  (D81+D  82p2)s«  ♦  (D71  ♦B72p2)s»         (3.12) 
♦  (D61  ♦D62P2)s6  ♦  (D51  *D52  p2)ss  +  (D41  +D42  p«)s* 
♦  (D  31*D^2  F2) s3  +  (D21*E22  P2)s2  ♦  D.,s  ♦  E^  =  0 


The  "migration"  of  the  system's  roots  as  the  rcll  rate 
increases  from  zero  is  shown  in  Figs.  3-3  through  3-12. 
Figure  3-4,  for  example,  shows  that  two  roots  of  the  system 
move  into  the  right-half  plane  at  a  rcll  rate  of  approxi- 
mately 12.2  radiar.s/second.  Since  the  linearized  character- 
istic equation  was  used  to  construct  the  root  loci  shown  in 
Figs.  3-3  through  3-12,  the  precise  value  cf  roll  rate  which 
causes  instability  cannot  be  determined  frcm  the  figures. 
Shat  the  figures  point  out  is  the  general  migration  cf  two 
system  roots  into  the  right-half  plane  in  all  phases  of 
operation  as  the  roll  rate  is  increased. 
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TABLE  3-1 


CHABACTEBISTIC  EQUATION  COEFFICIENTS 


1st  Stage 

2nd 

Stage 

3rd  Stage 

Beginning 

End 

Beginning 

End 

Beginning 

D9 

354.0 

354.0 

354.0 

354.0 

354.0 

°81 

39,879.0 

39,879.0 

39,879.0 

39,879.0 

39,879.0 

D82 

0.9755 

0.9755 

0.9401 

0.9401 

0.8869 

^L 

1.5809x10 

1.5809xl06 

1.5809xl06 

1.5809xl06 

1.5809xl06 

D72 

345.33 

345.33 

332.80 

332.80 

313.96 

°61 

3.0355xl07 

7 
3.0725x10 

3.0368xl07 

3.0751xl07 

3.0549x10 

D62 

38,902.0 

38,902.0 

37,490.0 

37,490.0 

35,368.7 

1   51 

3.1199xl08 

3.7781xl08 

3.1434xl08 

3.8236x10 

3.4650xl08 

1 

D52 

1.542xl06 

1.542xl06 

1.486xl06 

1.486xl06 

1.402x10 

D41 

1.7182xl09 

3.3458xl09 

1.7764xl09 

3.4584xl09 

2.571xl09 

°42 

7 
2.9483x10 

7 
2.9483x10 

2.8413x10 

2.8413xl07 

7 
2.6805x10 

D31 

4.8252xl09 

1.8385xl010 

5.310xl09 

1.9324xl010 

1.1929xl010 

°32 

2.815xl08 

2.815x10 

2.7128xl08 

2.7128xl08 

8 
2.5593x10    i 

D21 

7.311xl09 

8.2761xl010 

8.5265xl09 

8.1567xl010 

3.389xl010 

D22 

l.llxlO9 

l.llxlO9 

9 
1.0708x10 

1.0708x10 

1.0102x10 

Dl 

5.794xl09 

9.5278xl010 

7.0162xl09 

9.292xl010 

3.5410X1010 

°P__ 

1.931xl09 

3.176xl010 

2.3387xl09 

3.0972xl010 

1.1803xl010 

39 


i 


s 


o 


o 


o 

s 


PI 


s 

7 


* 


Irv^S) 


I 


fcw 


•IM.I   -1404       -1*0.0        -1004        -00.0         -00.0 


T f 

eo.o    0.0    to.o 


T 


T 


-404 


Fig-  3-3   Root  Migration  at  Beginning  of  First  Stage 

Figures  3-3,  3-5,  3-7,  3-9  and  3-11  show  how  the  roots 
of  the  characteristic  equation  of  the  celling  missile  move 
as  the  roll  rate  is  increased.  A  close-up  view  of  the 
migration  around  the  imaginary  axis  is  shown  in  Figs.  3-4, 
3-6,  3-8,  3-10  and  3-12. 
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Fig.  3-6   End  of  the  First  Stage  Boot  Migration  Close-Up 


43 


e 
e 

Zrr>(t) 



s 

© 

© 

$  i 

1   f?«ca; 

o 
© 

i 

*  \ 

M 
1 

o 

f- 

o 
1 

o 

§ 

— — — — 

-IMJ  -H04»    -IWLO    -100.0    -00.0    -004    -40.0    -MjO      0,0      tO.O 


Fig,  3-7   Root  Migration  at  Eeginning  of  Second  Stage 
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Fig.  3-9   Root  Migration  at  End  of  the  Second  Stage 
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Fig.  3-11   Root  Migration  at  Beginning  of  Third  Stage 
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B.   EFFECT  OF  BOLL  Oil  SISTEfi  PEBFOBHABCE 

As  shown,  the  introduction  of  roll  into  the  pitch-yaw 
controller  caases  a  migration  of  two  system  roots  into  the 
right-half  of  the  s-plane.  As  the  roots  approach  the  inagi- 
nary  axis  (as  the  roll  rate  is  increased)  the  system  begins 
to  show  more  transient  oscillation  due  to  the  decreased 
damping  ratio.  The  settling  time  of  the  system  is  also 
increased-  Since  the  system  contains  the  nonlinear  satura- 
tion effect  of  the  nozzle  actuator  and  is  of  such  high  cider 
(tenth),  not  much  more  than  this  qualitative  observation  was 
done  analytically. 

To  study  the  effect  cf  roll  on  system  performance,  a 
series  of  simulation  studies  was  conducted  with  the  satura- 
tion effect  present.  The  simulation  studies  were  performed 
by  subjecting  the  system  to  three  sets  of  step-input 
commands  and  constant  roll  rates.  The  first  study  consisted 
of  giving  the  system  a  step-pitch  input  and  a  zero-yaw  input 
and  observing  the  transient  response  at  six  roll  rates 
between  2.5  and  15  radians/second.  The  results  of  these 
simulations  are  plotted  in  Figs.  3-13  through  3-18.  The 
second  study  involved  subjecting  the  system  to  a  zero-pitch 
input  and  a  step-yaw  input  and  conducting  the  same  tests  as 
the  first  study.  The  results  of  the  second  study  are 
plotted  in  Figs.  3-19  through  3-24.  Finally,  the  system  was 
given  both  a  step-pitch  command  and  a  step-yaw  command  and 
the  tests  were  repeated.  The  responses  of  the  system  to  the 
conditions  of  the  third  test  are  shown  in  Figs.  3-25  through 
3-30.  Unlike  the  non-rolling  simulations  which  were  run  at 
separate  booster  stages  (Figs.  2-9  through  2-11),  the 
rolling  simulations  were  run  over  several  consecutive 
booster  stages.  This  is  because  the  settling  time  cf  the 
rolling  missile  system  was  often  greater  than  the  burn  time 
of  a  single  stage. 
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Fig.  3-13   Pitch-Step  JResponse  at  Boll  Eate  of  2.5  rad/sec 

figures  3-13  through  3-18  show  the  systei^s  pitch  and 
yaw  responses  to  a  pitch-step  input  and  a  zero-yaw  input  at 
the  specified  roll  rate.  Note  the  dramatically  increased 
overshoot,  oscillation  and  settling  tine  compared  to  the  the 
non-rolling  simulations,  Figs.  2-9  through  2-11-  figure 
3-18  shows  shows  the  unstable  response  obtained  at  a  roll 
rate  equal  of  15  rad/sec. 
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Fig.    3-14      Pitch-Step   Response   at    Boll    Sate    of    5.0   rad/sec 
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Fig.    3-15      Pitch-Step   Response    at    Eoll    Bate    of    7.5   rad/sec 
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Fig.    3-16      Pitch-Step    Response    at    Roll    Bate   of    10-0    rad/sec 
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Fig.    3-17      Pitch-Step    Response    at    Boll    Hate   of    12.5   rad/sec 
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Fig.  3-18   Pitch-Step  Eesponse  at  Eoli  Fate  of  15.0  rad/sec 
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Fig.  3-19   Yaw-Step  Response  at  Boll  Bate  of  2.5  rad/sec 

Figures  3-19  through  3-2<4  show  the  system's  pitch  and 
yaw  responses  to  a  yaw-step  input  and  a  pitch-zero  input  at 
the  specified  roll  rate.  In  these  simulations,  as  in  the 
previous  ones,  the  overshoot  and  settling  tine  are  increased 
above  that  observed  in  the  ncn-rolling  simulations  of  Figs. 
2-9  through  2-11.  Figure  3-2M  shows  the  unstable  response 
at  a  roll  rate  of  15  cad/sec. 
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Fig.    3-20      Yaw-Step    Response    at    Boll    Fate   of    5-0   rad/sec 
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Fig,    3-21      Yaw-Step   Response    at    Boll    Bate   of    7.5    rad/sec 
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Fig.    3-22      law-Step    Response   at   Roll    Bate   of    10.0    rad/sec 


60 


TIME  (SEC) 


Fig.    3-23      Yaw-Step    Besponse    at   Roll   Bate    of    12.5   rad/sec 
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Fig.  3-24   Yaw-Step  Besponse  at  Boll  Bate  of  15.Q  rad/sec 
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Fig.  3-25  Pitch-  8  Yaw-Step  Response  Foil  Bate  2.5  rad/sec 

Figures  3-25  through  3-30  show  the  system's  pitch  and 
yaw  responses  to  pitch-  and  yaw-step  inputs  at  the  specified 
roll  rate.  Under  these  conditions,  the  under-damped  oscil- 
lation observed  in  the  two  previous  test  is  even  nore 
pronounced. 
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Fig.    3-26   Pitch-    S    Yaw-Step    Response    Eoll    Bate    5.0   rad/sec 
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Fig.    3-27   Pitch-    6    Yaw-Step    Besponse    Boll    Bate   7.5    rad/sec 
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Fig.    3-28   Pitch-    6    Yaw-Step    Response    Roll    Fate    10.0   rad/sec 
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Fig.  3-29  Pitch-  £  Yaw-Step  Response  Roll  Bate  12.5  rad/sec 
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Fig.  3-30  Pitch-  6  Yaw-Step  Response  Boll  Bate  15.0  rad/sec 
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The  simulation  studies  showed  that  the  systen  begins  to 
reach  instability  at  a  roll  rate  of  approximately  12 
radians/second-  However,  at  roll  rates  as  low  as  2.5 
radians/second  the  system's  performance  has  been  severely 
degraded  from  the  non-rolling  case  {figs.  2-9  through  2-11). 
In  all  cases,  the  maximum  peaX  overshoot  is  increased  by  a 
factor  of  at  least  five  and  the  settling  time  by  a  factor  of 
at  least  eight.  The  system's  performance  is  approximately 
that  of  an  under-damped  oscillator.  This  response  is  caused 
by  tie  two  roots  which  have  moved  close  to  the  imaginary 
axis  as  the  roll  rate  was  increased  and  are  dominating  the 
system's  transient  response.  In  the  next  section,  a  compen- 
sation scheme  is  proposed  which  increases  the  damping  of  the 
system  and  partially  cancels  the  cross-coupling  acceleration 
terms  introduced  ty  roll. 
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IV.  COHPEBjSAlION  SCH|fl| 

In  the  preceding  chapter  it  was  noted  that  rolling  the 
missile  had  the  effect  of  introducing  cross-coupling  terms 
into  the  pitch  and  yaw  angular  acceleration  eguations  (Eqs. 
(3.6  and  (3.7)).  In  order  to  completely  restore  the 
system's  performance  to  its  non-rolling  condition,  it  is 
necessary  to  eliminate  or  cancel  out  the  cross-coupling 
terms.   Equation  (3.10)  can  te  rearranged  into  the  form: 


s* 

Qs' 

[9  1  = 

fG<9c    * 

re 

-    9)" 

es 

s*_ 

w 

H(fc    - 

-    KrfS 

-    Y ) 

(4.1) 


The  left  side  of  Eg.  (U.  1)  contains  the  "dynamics"  of  the 
system.  In  other  words,  it  gives  the  acceleration  and 
velocity  terms  which  determine  the  angular  acceleration  in 
the  pitch  and  yaw  channels.  Ihe  only  way  to  eliminate  the 
cross-coupling  terms  directly  would  be  to  make  changes  to 
the  left  side  of  Eg.  (4.1)  which  would  cancel  the  terms,  Qs 
and  Es6.  This  would  involve  dramatic  charges  to  the  missile 
body,  thrust  vector  control  system  and  flight  performance 
and,  therefore,  will  not  be  attempted. 

The  right  side  of  Eg.  (4.1)  constitutes  the  "elec- 
tronics" of  the  system.  It  contains  a  sum  of  those  terms 
which  are  added  at  the  feedback  junction  then  passed  through 
the  nozzle  actuator.  This  side  of  the  eguation  can  be 
changed  easily  by  altering  the  electronic  network  that  makes 
up  the  feedback  lcop. 
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A.   CBOSS-CHASNEI  FEEDBACK 

Assume  that  the  pitch  rate  and  yaw  rate  can  be  accu- 
rately measured  or  estimated  at  any  tine  during  flight. 
Also,  assume  that  the  terms,  6;  H;  C  -and  E,  all  functions  of 
system  parameters,  can  be  measured  or  estimated  accurately. 
Then,  addition  of  the  feedback  terms  cf  pitch-rate-to-yaw- 
input  and  yaw-rate-to-pitch-input  permit  cancellation  of  the 
cross-coupling  terms.  The  block  diagram  realization  of  the 
proposed  change  is  shown  in  Fig.  4-1.  Equation  (4.2)  gives 
the  compensated  system's  transformed  differential  equation 
in  matrix  form. 
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Equation  (4.2)  can  be  reduced  to  Eg  (4.3)  by  canceling 
the  cross-coupling  terms,  Qs  and  PsC,  from  both  sides  of 
the  eguation.  Comparison  of  Eq.  (4.3)  with  Eg.  (3.5)  shows 
that  the  transformed  differential  equation  of  the  rolling 
system  has  been  reduced  to  its  original  non-rolling  form. 
Note  that  all  that  was  done  in  this  compensation  scheme  was 
to  subtract  the  cross-coupling  term  from  the  command  input 
to  the  nozzle  actuator  so  than  when  it  is  added  back  on 
downstream  the  result  will  be  no  effective  cross-channel 
terms  into  either  ccntrol  path. 
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Fig.  4-1   Cross-Channel  Feedback  Compensation  Scheie 
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Unfortunately,  the  problem  is  not  as  simple  as  Eg.  (4.2) 
implies-  The  main  problem  is  that  the  amount  of  control 
which  can  be  used  to  force  the  system^  response  is  limited 
by  the  maximum  deflection  angle  of  the  eihaust  noz2le.  This 
nonlinear  saturation  of  the  actuator  is  not  present  in  Eq. 
(4.2) .   Figure  4-2  summarizes  the  problem. 
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Fig.  4-2   Plant  Dynamics  and  Control  Input  Limitations 

Figure  4-2  shows  that  no  matter  what  feedback  scheme  is 
employed  to  generate  the  pitch  and  yaw  angular  acceleration 
control  signals,  they  are  limited  and,  therefore,   so  is  the 
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region  of  operation  over  which  the  performance  of  the  system 
can  be  controlled  and  the  degree  of  control  possible. 

Figure  4-3  shows  the  cross-coupling  terms  and  the 
attempted  cancellation  of  them  by  the  cross-channel  feedback 
terms  under  the  compensation  scheme  just  proposed.  The 
simulation  shown  was  performed  during  the  first  stage  of 
operation  with  the  system  subject  to  a  pitch-step  command 
and  a  zero-yaw  command.  Note  how  the  magnitude  of  the 
cross-channel  feedback  term  is  effectively  saturated  by  the 
nozzle  actuator  and  cannot  folio*  the  cross-coupling  term 
outside  a  small  region  of  operation.  Ihe  simulation  was 
made  by  subtracting  the  error  tera  generated  by  the  differ- 
ence in  command  and  feedback  so  that  only  the  cross-channel 
feedback  term  is  fed-through  the  nozzle  actuator. 

B.   NONLINEAB  FEEDBACK 

Examination  of  Figs.  2-5,  2-9,  2-10  and  2-11  shows  that 
the  time  the  nozzle  actuator  takes  to  reach  maximum  deflec- 
tion when  given  a  large  input  command  (>>  maximum  deflec- 
tion) is  insignificant  (0.05  sec)  when  compared  to  the 
settling  time  of  the  entire  pitch-yaw  contrcl  system.  Thus, 
the  actuator  is  functioning  almost  like  a  switch  in  all 
regions  of  operation  except  near  the  coamand  input  (steady 
state).  For  example.  Fig.  4-4  shows  the  nozzle  deflection 
in  the  pitch  plane  during  a  simulation  when  the  systea  was 
subject  to  a  pitch-step  input  and  a  roll  rate  of  2.5 
radians/second.  Ihe  conclusion  to  be  drawn  from  Fig.  4-4  is 
that  since  the  aaplitude  of  the  effective  control  input  to 
the  plant  cannot  be  increased  above  saturation,  the  only 
alternative  is  to  "throw  the  switch"  a  little  sooner. 

The  cross-channel  feedback  compensation  scheae  just 
proposed  will  cause  a  reduction  in  the  effective  servo 
command  and  result  in  earlier  removal  of  the  aaxiaum  nozzle 
deflection  (turning  off  the  switch).  However,  this  approach 
can  at  best  only  follow  the  cross-coupling  teras  and  can  not 
"anticipate"  the   need  to  switch   the  deflection   angle  from 
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positive  maximum  to  negative  maximum.  Another  method  to 
accomplish  the  objective  of  switching  the  nozzle  earlier  is 
to  vary  the  velocity  feedback  gains,  K  re  and  Kro/  . 
Specifically,  it  is  desired  that  the  velocity  feedback  gains 
fce  low  during  the  initial  phases  of  operation  when  the  error 
between  command  and  position  is  large  and  increase  as  the 
system  approaches  the  commanded  value.  This  will  result  in 
the  same  basic  system  operation  that  was  present  before 
compensation  {dominant  complex-pair  oscillation)  but  with 
low  damping  when  the  error  is  high  and  higher  damping  as  the 
system  approaches  steady  state. 

The  nonlinearity  of  the  system  prevents  the  use  of  most 
analytical  tools  which  might  be  used  to  determine  algorithms 
for  K  r0  and  Kry  .  Instead,  a  trial  and  error  approach  was 
used  based  on  the  observations  of  the  system^  performance 
during  earlier  simulations  and  simulation  studies  performed 
specifically  to  determine  appropriate  parameters. 

Assume  that  the  algorithm  for  determining  the  velocity 
feedback  gain  consists  of  a  constant  which  is  reduced  by  an 
amount  proportional  to  the  error  between  command  and  posi- 
tion. This  will  result  in  low  damping  if  the  error  is  high 
and  increased  damping  as  the  error  decreases.  For  example, 
if  V  and  H  are  constants,  the  pitch-channel  velocity  feed- 
back gain  can  be  determined  from: 


Kre  =  V  -  W  9c  '  9 


(4.4) 


Furthermore,  the  simulations  shown  in  Figs.  3-13  through 
3-30  indicated  that  the  system  becomes  more  underdamped 
(roots  migrate  toward  the  imaginary  axis)  as  the  roll  rate 
is  increased.  So,  the  velocity  feedback  gains  were  made 
proportional  to  the  coll  rate.  Normalizing  the  algorithm  so 
that  the  relationship  does  not  become  zero  if  tbe  roll  rate 
is  zero  yields  Egs.  (4.5)  and  (4.6). 
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K         =     (V    -    M    0       -    9    )  {1    ♦    p) 

re  u  i 


C-5) 


=    (V   -    it     -  Y    )  (1    *    F) 


(4.6) 
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Fig.  4-3   Saturation  of  Cross-Channel  Feedback 
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Fig-  4-4   No2zle  Deflection  Angle  Siiulation 
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Figure  4-5  shows  sample  results  of  simulation  studies 
which  were  conducted  for  various  values  of  V  and  H  and  ever 
various  stages  of  operation  and  roll  rates.  The  simulations 
showed  that  acceptable  performance  was  obtained  for  V  =  4 
and  N  =  3. 

Figures  4-6  through  4-14  show  simulation  results  which 
were  made  by  using  the  cross-channel  feedback  and  variable 
velocity  feedback  gain  coapensation  schemes  together.  For 
comparison,  the  uncompensated  responses  are  also  shown.  The 
compensated  simulation  responses  were  obtained  under  the 
same  operating  conditions  as  those  shown  in  Figs.  3-13 
through  3-30. 


V,W)  -  (2,1) 
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Fig.  4-5   Sample  Algorithn  Eesponse  Test  Simulation 
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Fig.    4-6      Pitch-Step  Besponse    Boll    Bate    of    2-5    cad/sec 
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Fig.  4-7   Pitch-Step  Response  Roll  Bate  of  5.0  rad/s€c 
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Fig-    <4-8      Pitch-Step   Response    Boll    Bate    of    7.5    rad/sec 
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Fig.    4-9      Yaw-Step   Response    Boll    Pate    of    2-5    cad/sec 
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Fig.    4-10      Yaw-Step    Response    Roll    Rate   of    5.0   rad/sec 
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Fig.    4-11    Yaw-St€p   Response    Foil    Rate    of    7.5   rad/sec 
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Fig.    4-12   Pitch-    S    law-Step    Besponse    Boll    Bate    2.5   rad/sec 
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Fig.    4-13    Pitch-    6    Yaw-Step    Response    Boll    Bate    5-0    rad/sec 
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Fig  4-14  Pitch-  6  Yaw-Step  Besponse  Bell  Bate  7.5  rad/sec 

The  siaulatioDS  indicate  that  substantial  performance  resto- 
ration is  possible  with  the  proposed  compensation  schenes. 
However,  the  restoration  is  not  complete  and  is  still  much 
worse,  in  terms  of  overshoot  and  settling  time,  than  the 
non-rolling  missile. 
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V.  CONCLUSIONS 

The  compensation  scheme  proposed  partially  restores  the 
pre-roll  performance  of  the  missile-  Table  5-1  summarizes 
the  level  of  performance  degradation  caused  by  roll  and  the 
restoration  which  was  obtained  with  this  compensation  scheme 
under  pitch-  and  yaw-step  inputs  and  a  roll  rate  of  2-5 
rad/sec.  As  indicated,  the  degradation  is  significant  even 
at  this  relatively  low  (and  stable)  roll  rate.  The  compen- 
sated system's  performance  is  still  far  worse  than  that  of 
the  non-rolling  system- 

The  compensation  scheme  relies  primarily  on  the  nonli- 
near damping  effect  of  the  variable  velocity  feedback  gains 
which  "anticipate"  the  approach  of  the  commanded  position. 
Since  the  cross-coupling  terns  introduced  by  roll  are  not 
directly  cancelled  by  the  proposed  modification,  the  region 
of  stable  operation  and  the  degree  to  which  the  system's 
performance  can  be  restored  is  directly  controlled  by  the 
roll  rate-  Furthermore,  the  saturation  caused  by  the  nozzle 
deflection  will  prevent  any  compensation  scheme,  which  only 
modifies  the  control  input,  from  fully  restoring  the 
system's  performance-  To  completely  restore  the  system  to 
its  pre-roll  performance  at  any  roll  rate,  additional 
control  inputs  to  the  system  plant  dynamics  (such  as,  side 
thrusters  or  additional  vectoring  of  the  main  rocket  motor) 
are  necessary- 
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1ABLE  5-1 


PERFORMANCE  COMPARISON 


Non-Rolling 
System 

_._,  . 1 

Uncompensated    Compensated 
Foiling  System     System 
(p  =  2.5  cad/sec) 

Settling  Tine 

9  sec 

140  sec 

80  sec 

Peak  Overshoot 

. .  _ 

8. 571 

170* 

70* 

__.,   .„,.,.,] 

89 


APPENDIX  A 
NON-BOLLING  HISSILE  SIHULA1ICN 

The  DSL/360  program  used  to  simulate  the  non-rolling 
missile  is  shown  below.  The  program  consists  of  three 
stages,  initial,  derivative  and  terminal.  The  initial 
segment  is  used  to  define  constants  and  calculate  functions 
of  the  constants  which  will  net  change  throughout  the  simu- 
lation run.  Ihe  derivative  section  is  where  the  actual 
integrations  reguired  to  perform  the  simulations  are 
conducted.  This  program  is  constructed  to  follow  the  block 
diagram  of  Fig.  2-1  directly.  Note  that  the  nozzle  actuator 
transfer  function  was  broken  down  into  a  three-step  integra- 
tion. Ihe  terminal  stage  determines  the  end  condition  of 
the  simulation  run  and  the  output  data.  This  data  was 
formed  into  a  data  file  and  a  packaged  graphics  program 
(DISSPLA)  was  used  to  construct  the  figures  shown  in  the 
text.  The  parameter  values  listed  in  this  example  program 
are  those  during  stage  one. 


*  INITIALIZE  CONSTANTS 

CONST  TST  =  7496000, L=363.  J  1=110*48. 8  .KTH=1.0,KBIH=  1.75, 
WA=150,WB=15,PB=0.9,BAIE  =  930.0,M=1891'40 

*  COMPUTE  SET  VALUES 
INITIAL 

CNE=WA*  (WB**2) 
IWO=HA  ♦  2*PB*WE 
IHBEE=2*PB*SB*WA  ♦  WB**2 
* 

*  BEGIN  SIHULATICN 
DEBIVATIVE 

HASS=M1-BATE*TIME 
PCCM=STEP(0) 

*  EBBOB  SIGNAL 

E=ECOM-P 

EI  =  E-  <KBTH*E1C0T) 

CM=KTH*EI 

*  SEBVO  TBANSFEB  FUNCTION 

SVO  =  INTGBL  IQ. SV1DOI1 

SV1D0T  =  INTGBI  (0,SV2DCT) 

SV2DOT=INTGBL  (0,SVJDGT) 

SV3DOT=ONE*  (CH-SVO)  -  TWO+SV2DCT  -  THBE£*SV1D0T 
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*  LIMIT  OF  NOZZLE  DEFLECTION 

DEJL=LIMIT <-0.0  873#0.08  73,SVO) 

*  MOMENT  OF  IREBTIA 

E=J1*MASS 
Z=TST*L/B 

*  ANGULAB  ACCELEBATION 

P2DOT=Z*DEFI 

*  ANGULAB  VELOCITI 

P1D0T=INTGRL (0,P2DOT) 
P=INTGBL  (0,ElDOI) 

*  OUTPUT  AND  CONTBOL  STATEMENTS 
IEBMINAL 

CONTBL  FINTIfl=15.0,  CELT  =  0.001 

SAVE   0.015,P 

PRINT  0.015,P 

END 

STOP 
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BOLLIHG  HISSILE  SI8ULATICN 

Unlike  the  non-rolling  system  simulations,  the  rolling 
missile  simulations  were  required  to  run  over  several 
consecutive  booster  stages.  This  required  the  addition  of  a 
•dynamic'  segment  of  simulation.  The  'IF*  statements,  shown 
in  the  example  sinulation  below,  are  used  to  switch  from  one 
set  of  parameters  to  another  as  one  stage  burns  off  and 
another  starts.  The  switching  is  controlled  by  the  elapsed 
time  of  the  simulation  run  compared  to  the  turn  time  of  each 
stage. 


*  SET  CONSTANTS 

CONST  KTB=1.0,KBTH=1. 5,WA=150, WB=  1  5, PB=0. 9# BCLL=7. 5 

*  INITIALIZE  SET  VALUES 
INITIAL 

A  =  HA*|HE**2) 

B  =  HA  ♦  2*FB*WE 

C  =  2*PE*WB*WA  +  WB**2 

E  ES ET  =  0  0 
♦INITIALIZE  FIBST  STAGE  VALUES 

TST  =  7<496000 

L  =  363 

J1=11048.8 

J2=136.  1 

EAIE=930.0 

Ml=189140 
* 

*  DECIDE  WHICH  STAGE  OF  FLIGHT 
DYNAHIC 

*  ENTEE  STAGE  THO 

IF  (TIME. LT. 150)  GO  IC  10 
TST=1125000 
L  =  230 
J1=4476.4 
RATE=81.49 
M1=40340 
RESET=150 

*  ENTEB  STAGE  THRIE 

IF  (TIME. LT. 509)  GO  TC  10 
TST=224900 
L=100 
J1=858.3 
J2=50.0 
EATE=14.75 
M1=8140 
BESET=509 
10     CONTINUE 


*  SIMULATION  BUN 
DEBIVAIIVE 

SET=TIME-BESET 

MASS  =  M1  -  <BATE*SET) 

PCCM  ■  STEP  (o; 

YCGM  -  STEP 

*  EBBOB  SIGNALS  POSITION  FEEDBACK 
EP  =  PCOM-P 
EY  =  YCGM-Y 

*  EBBOB  SIGNALS  VELOCITY  FEEDBACK 
EIP  =  EP-(KBIH*P1D0I: 
EIY  = 
CMP  = 

CHY  =  K1B*EIY 

*  SERVO  TBANSFEB  FUNCTION  IN  PITCH  PATH 


'J°l 

»CSI' 
I 

I 

t/  ELOC  j.  ix  re 
EP-(KBTH*P1D0I) 
EY- JKBTH*Y1DOT) 
KTH*EIP 


SVP  =  INTGBL(0,SVPDOT) 
SVPDOT  ■  INTGBL  I0,SVP2DT) 
SVP2DT  =  INTGBL  (0#SVP3DT) 


SVP3DT  =  A*<CMP-stfP)  -  E*SVP2DT  -  C*SVPDOT 

*  PITCH  DEFLECTION 

PDEFL  -    LIMIT  1-0.0873- 0.0873. SVP) 

*  SEBVO  TBANSFEB  FUNCTION  IN  YAH  PATH 

SVY  =  INTGBl(O.SVYDOT) 

SVYDOT  =  INTGBL <0,SVY2DT) 

SVY2DT  =  INTGBL  |0,SVY3DT) 

SVY3DT  =  A*(CMY-SVY)  -  E*SVY2DT  -  C*SVYDOT 

*  YAH  DEFLECTION 

YDEFL  =  LIMIT J-0. 0873, 0.0873, SVY) 

JY  =  J1*HASS 

JZ  =  JY 

JX  =  J2*MASS 

ZP  -  TST*L/JY 

ZY  =  TST*L/JZ 

*  CBOSS  COUPLE  MULIIPLIEBS 

ALPHA  =  (JZ-JX) *BOLL/JY 
EETA  =  <JX-JY)*BOLL/JZ 

*  CBOSS  COUPLE  TEBH  YAH  TO  PITCH 

CCYP  =  ALPHA*Y1D0T 

*  CBOSS  COUPLE  TEBH  PITCH  TO  YAH 

CCPY  =  BETA*P1DOT 

*  PITCH  ANGULAR  ACCELERATION 

P2DOT  =  (ZP*PDEFL)  ♦  CCYP 

*  YAH  ANGULAB  ACCELEBATION 

Y2DOT  =  (ZY*YD£FL)  ♦  CCPY 

*  PITCH  ANGULAB  VELOCITY 

P1DOT  =  INTGBL (0.P2DOT) 
P  =  INTGBL  (0.P1DOT) 

*  YAH  ANGULAB  VELOCITY 

Y1DOT  =  INTGBLI0.Y2DOT) 
Y  =  INTGBL (0,I1DCT) 
* 

*  OUTPUT  AND  CONIBOL 
TERMINAL 

CONTBL  FINTIH=500.0,  DELT=0.01 

PBINT   0.50,P,Y 

END 

STOP 
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